Liouvillian Approach to the Integer Quantum Hall Effect Transition 
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We present a novel approach to the localization-delocalization transition in the integer quantum 
Hall effect. The Hamiltonian projected onto the lowest Landau level can be written in terms of the 
projected density operators alone. This and the closed set of commutation relations between the 
projected densities leads to simple equations for the time evolution of the density operators. These 
equations can be used to map the problem of calculating the disorder averaged and energetically 
unconstrained density-density correlation function to the problem of calculating the one-particle 
density of states of a dynamical system with a novel action. At the self-consistent mean-field level, 
this approach yields normal diffusion and a finite longitudinal conductivity. While we have not 
been able to go beyond the saddle point approximation analytically, we show numerically that the 
critical localization exponent can be extracted from the energetically integrated correlation function 
yielding v — 2.33 ± 0.05 in excellent agreement with previous finite-size scaling studies. 
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I. INTRODUCTION 

The metal-insulator transition in the integer quantum 
Hall effect (IQHE) is a reentrant zero temperature quan- 
tum phase transition in which the sample goes from an 
insulating phase with longitudinal conductivity a xx = 
to another insulating phase by crossing a conducting crit- 
ical point (cr xx 7^ 0) as the magnetic field is varied. The 
critical point occurs between the plateaus of the Hall con- 
ductivity cr xy and corresponds to the instance when the 
Fermi energy is at a critical energy located in thejniddle 
of one of the disorder broadened Landau levels.ErH 

In general, the disorder induced metal-insulator tran- 
sition is a transition in the nature of the states (whether 
they are localized or delocalized) at the Fermi energy 
and it does not manifest itself in the density of states 
which remains smooth across the mobility edge. Accord- 
ing to the one-parameter theory of scaling, the states of 
a two-dimensional noninteracting electron gas are all lo- 
calized in the presence of arbitrary weak disorder .0 In the 
IQHE however, the presence of the strong magnetic field 
pointing perpendicular to the plane drastically changes 
the nature of the states near the middle of the Landau 
bands. In the noninteracting picture of the IQHE these 
states are characterized by a localization length 



E- El 



E 



(1) 



which determines the extent to which the eigenstates of 
energy E are delocalized. Here £o denotes a characteristic 
length scale of the system, e.g. the magnetic length t 
(see below) and Eq a characteristic energy scale, e.g. the 
bandwidth or disorder strength. The critical energy E l c 
is located in the middle of the i-th Landau band and, 



in an infinite size system, it is the only energy at which 
the one-particle eigenstates are delocalized within this 
Landau band. As the Fermi energy (or magnetic field) 
is varied, the conductivity c xx will change according to 
the nature of the states at that energy and sharp peaks 
in the longitudinal conductivity will be observed. 

When studying the IQHE the interaction between the 
electrons is usually ignored and only the disorder is con- 
sidered to be responsible for the localization of the sin- 
gle particle states. This assumption must be checked by 
comparing the prediction^, of the noninteracting theory 
to experimental resultsatm and the outcome of-rwmer- 
ical calculations which include the interactions.B~l3 The 
universal localization exponent v — 2.34 ± 0.04 numer- 
ically obtained within a noninteracting theoryB~t3 is in 
excellent agreement with experimental measurements of 
i/J^'Cru but it remains a mystery why the strong interac- 
tions, which do afied-the dynamical exponent z, does not 
seem to affect zAtnDO Here we adopt the noninteracting 
picture. We furthermore assume a strong magnetic field 
and a Zeeman splitting, which is much larger than the 
width of each disorder broadened Landau level. We can 
then focus on the transition within the lowest Landau 
level (LLL) and neglect the spin degree of freedom of the 
electrons. 

It has been shown numerically that for a finite sys- 
tem delocalized one-particle wave functions near E c show 
multifractal properties characterized by a set of general- 
ized fractal dimensions Z5 q .EHl3't3 Also, dynamical studies 
have shown anomalous slow diffusion of-wave packets con- 
structed from these multifractal states .til Diffusion can be 
studied using the spectral function of the disorder-aver- 
aged retarded density-density correlation function.ll3 For 
the problem considered hcrp-the spectral function is given 
by (after dividing by nhujjEl 
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S(r,w;E) = (ly^diE -Jicj/2- Ei)S(E + ncj/2- Ej) 



1,3 



x^(0)V*(r)^-(r)V*(0))). 



(2) 



Here the ipi( r ) denote one-particle eigenfunctions and 
Ei the respective eigenenergies for an electron of a two- 
dimensional spinless electron gas which is subject to a 
perpendicular magnetic field and a disorder potential. 
((. . .)) indicates the ensemble average over the disorder. 
After taking the disorder average, translational invari- 
ance is restored and S only depends on the distance 
r = |r| from the origin of the plane. Assuming that the 
eigenstates which contribute in Eq. (0) for E E c are of 
multifractal piia|racter, it has been argued that S decays 
algebraicallyolla 



S(r,w-+ 0;E^ E c ) 



(3) 



for £q ^ r ^ The anomalous diffusion exponent r\ 

is related to the generalized fractal dimension via D2 = 
2 - rjm 

Assuming a generalized nonlocal (in time and space) 
relation between the current and the gradient of the den- 
sity and using the continuity equation, the spectral func- 
tion in momentum space S{q,u>;E) at small q = |q| and 
u>, can be rewritten in terms of— a, generalized diffusion 
"coefficient" D(q,u>) for E « EJB 



S(q,u;E) 



P(E) 



hq 2 D(q,to) 



[hu] 2 + [hq 2 D(q,l 



(4) 



where p{E) is the density of states per unit area. In 
the limit of u>, q — > and for large enough system sizes, 
D(q,u>) is only a function of qL u , where£3 



= [p(E c )hu>] 



-1/2 



(5) 



Through numerical diagonalization and using Eq. (Q) 
Chalker and DanielltHl have shown that D(q,u>) ap- 
proaches a constant Dq for small qL u . The precise value 
of -Do is important since the longitudinal conductivity 
at the critical point is given by the Einstein relation 
<T XX = e 2 p(E c )Dq and is expected to be universal. The 
qL u —> limit of D(q, u>) has later been reinvestigated in 
an extended numerical study.EH For qi^r^ 1, but still in 
the limit of q, uj — > 0, D(q, us) decays ast2l 



D(q,u) cx D (qL^-v . 



(6) 



For the anomalous diffusion exponent r\ Chalker and 
Daniell obtain the numerical value r\ = 0.38 ± 0.04, indi- 
cating that the delocalized states near the critical energy 
indeed have multifractal properties. This value fofpf* has 
later been confirmed in other numerical studies. L-HtLS For 
energies E away from E c , S(q, u>; E) vanishes in the small 
q and uj limit independent of the order in which the limits 
are taken due to exponential localization of the states. 



Most of the progress in the theoretical understand- 
ing of the localization-delocalization transition coin 
ered here has been through numerical calculationsflu 
Although a field theory has beea proposed some time 
ago by Pruisken and co-workers,t£l up to now no quanti- 
tative results such as the critical exponents of the tran- 
sition have been obtained within this description. More 
recent studiesoEj have introduced alternative field the- 
ories. Within the framework of these theories it might in 
the future be possible to analytically determine critical 
exponents as has been recently successfully achieved for 
the SU(2) version of the network modelB3 In this paper 
we present a novel approach to the transition which may 
prove more tractable. Although thus far we have not 
been able to analytically calculate the spectral function 
Eq. (||) beyond the self-consistent Born approximation, 
we have numerically verified the possibility of obtaining 
the critical exponent v using this approach. 

We start by defining the density correlation function 
at zero temperature as 



i6(t) 
NM 2 



«Tr{p q (t)p_ q (0)«5(£ - H)})) , (7) 



with the one-particle Hamiltonian H = H + Ho- Here 
Hq denotes the kinetic energy of a spinless electron mov- 
ing in the plane in the presence of a perpendicular mag- 
netic field and Hd is the potential energy for a fixed real- 
ization of the disorder potential V(r). I is the magnetic 
length given by £ 2 — hc/(eB), where B is the strength 
of the magnetic field, and N = L 2 /(2irl 2 ) is the num- 
ber of states in the LLL. We consider a square sample of 
area L 2 . By projecting the one-particle density operator 
Pq = exp (— iq ■ r) onto the LLL, denoting the projected 
density by p q (see Sec. II), and taking the one-particle 
trace Tr over the states in the LLL, we restrict our con- 
siderations to the transition in the LLL. It will turn out 
that the equation of motion for the density operators re- 
stricted to the LLL can be solved formally in this case. 
In the small lu limit we have 



S{q,u>;E) = -^ImIl(q,Lj;E). 



(8) 



Instead of dealing with Yi(q,t;E w E c ) directly, we will 
integrate H(q, t; E) over all energies E and focus our at- 
tention on 

U(q,t) = - l -&((Tr{p- q (i)p- q (0)})) . (9) 

Since the localization length only diverges at E c , the en- 
ergetically unconstrained diffusion problem considered by 
investigating Imn(q,w) still contains useful information 
about critical exponents. For instance, let us suppose 
that at time t — we create a wave packet localized at 
the origin constructed from all the states of the system 
(localized as well as delocalized states). For large t only 
"delocalized" states with £(E)/r > 1 can contribute to 



2 



the probability amplitude of the wave packet at a dis- 
tance r/£o 3> 1 far away from the origin. This implies 
that in the limit of small q£o only states with q£(E) > 1 
and thus \E — E c \< E (q£ Q ) 1 / v contribute to the right 
hand side of Eq. ( |) . Hence for q£o — * only a fraction 
~ (qCo) 1 ^ of the states in the LLL contributes and we 
expect from Eq. (0), that for small qL^ 



ntui 2 lmll(q 7 uj) cx (q£o) 



xlv (D q 2 M 
1 + (D q 2 H 2 



(10) 



where the diffusion parameter is a constant Dq. The 
above argument, which we confirm numerically in Sec. 
V, gives a strong indication that some useful information 
about the quantum phase transition can be extracted 
from H(q,uj). We again emphasize that this is so be- 
cause the derealization only occurs at a single critical 
energy, a characteristic unique to the IQHE where the 
extended states have zero measure in the energy spec- 
trum. We also point out the importance of the order 
of limits in obtaining Eq. ([To|). The limit of q,uj — > 
is taken by having q approach zero faster than u> so as 
to obtain a finite diffusion constant. In contrast to the 
usual approach,E3 in which information about the anoma- 
lous diffusion exponent r\ is extracted from the spectral 
function S(q,uj; E rs E c ) we will be able to extract infor- 
mation about the localization exponent v using the same 
spectral function but integrated over all energies E. 

We will show that H(q, uj), which is an inherent 
fermionic disorder averaged two-particle correlation func- 
tion, can be re-expressed as the single particle correlation 
function of an interacting (after the disorder average has 
been performed) dynamical system with an unusual ac- 
tion. Therefore, in order to extract the dynamical behav- 
ior of the original problem, one simply has to study the 
disorder averaged density of states of this new action. 

The rest of this paper is organized as follows. In Sec. II 
we introduce the model and mapping of the problem to 
the new "Hamiltonian" . In Sec. Ill we calculate H(q, cu) 
within the self-consistent Born approximation. It dis- 
plays normal diffusion at this level of approximation. In 
Sec. IV we introduce the field theoretical approach to the 
disorder averaging. In Sec. V we will demonstrate numer- 
ically the validity of the scaling hypothesis stated in Eq. 
(|l0|), and finally in Sec. VI we present our conclusions. 



II. MODEL AND MAPPING 

We consider the two-dimensional spinless electron gas 
lying in the x-y plane which is subject to a perpendicular 
magnetic field B = Bi and an external potential V(r). z 
denotes the unit vector in the z direction. In the symmet 
ric gauge the vector potential is given by A : 
and the one-particle Hamiltonian reads 



irxB 



H 



H - 
1 
2^ 



P + 



e -A 



V(r) 



(11) 



We restrict our investigations to the LLL and thus 
project the Hamiltonian onto the states in the LLL. The 
kinetic energy of all the LLL states is the same and af- 
ter projecting leads to a constant which we will neglect 
in what follows. Writing the potential energy in Fourier 
space the Hamiltonian simplifies to 



H = ^v(-q)p q 



(12) 



where w(q) is the Fourier transform of the disorder po- 
tential. The projected density operator is given by 



Pq 



'q ) 



(13) 



with r q being the unitary magnetic translation operator 
which translates the electron a distance (q x z)£ 2 . The 
formalism needed to project the density operator p q = 
e -jq r on to the LLL was developed elsewhere.cJ 

The magnetic translation operators have the following 
special property: 



T a T„ 



exp 



-qAp 



'q+p i 



(14) 



where q A p = (q x p) ■ z. Hence their commutation 
relation defines a closed Lie algebra: 



[rq, t p ] = 2isin ( — q Ap ) r q+p 



Also we have 



Tr{r q } = iV5 q , . 



(15) 



(16) 



The latter can be proved by noting that the left hand 
side is proportional to the one-particle trace of p q . Since 
the trace is taken over states in the LLL, the projection 
is unnecessary and we have 



Tr{p q } = Tr{e-^ r } 



(17) 



which vanishes unless q = 0. 

If there are N states in the Hilbert space, there are 
A^ 2 independent operators on the space. However there 
are exactly iV 2 different wave vectors on the torus, so 
the set of operators p q is "complete"; it spans the set 
of all operators. The Hamiltonian can be expressed in 
terms of the p q and the Heisenberg equation of motion 
of the p q is closed. This allows us to define the quantum 
"Liouvillian" matrix by 



■ ^qq'Tq'(t) 



From the simple commutation properties Eq. (|lf 
r q it readily follows that 



(18) 



of the 



2i 

T 



«(q- q')e" 



-q| 2 



q'Aq). (19) 
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Using the Liouvillian matrix we can immediately write 
down the formal solution of the equation of motion Eq. 
(H) for r q (t) 



7q(*) 



q' 



Mo). 



(20) 



This leads to a simple expression for the density-density 
correlation function defined in Eq. (0) 



U(q,t) 



At) - 

-i — — — p 

hi 2 



-iCt\ 



' qq 



(21) 



We can define an iV 2 element operator "superspace" and 
view C as the "Hamiltonian" . From this point of view, 
finding ImII(g, u>) is the same as finding the one-particle 
density of states for a system with Hamiltonian C: 



n(g,w) 



p CO 

' le2 <> 2 / rfte^+ i5 )'(((q| 
Jo 



M 2 



e 2 



-iCt 



M 2 



1 



lu + iS — C 



q>» 



(22) 



where we have introduced states |q) with (q| C |q') = 
£ qq / and S is an infinitesimal small positive number. 

This remarkable formula is our central result. Let us 
now try to understand its import. In a crude sense it 
represents a kind of bosonization of the problem. Ordi- 
narily in an interacting many-body system the equations 
of motion for the density are not closed but rather in- 
volve a hierarchy of additional operators. However for 
the special case of one-dimension and a linear dispersion 
relation (the Tomonaga-Luttinger model) the equations 
of motion ape. closed and the density fluctuations become 
free bosonsE-3 even though the underlying particles are 
interacting. In the present problem (without electron- 
electron interactions) the equations for the density op- 
erators close after projection onto a single Landau level 
(which for simplicity we have taken to be the lowest). 
This has several advantages. First we do not have to 
work separately with retarded and advanced one-particle 
Green's functions and their products. Secondly we note 
that there are no problems with gauge invariance and 
conserving approximations. This is because the Liouvil- 
lian matrix elements £ qq ' vanish if either q or q' vanish. 
Thus the total charge in the system is automatically con- 
served. Finally this representation allows us to establish 
a hierarchy of length and time scales which should be 
suitable for renormalization group (RG) analysis. Be- 
cause the kinetic energy has been quenched, high mo- 
mentum of a particle is not associated with high energy. 
Since the Liouvillian vanishes at small wavevectors, it 
naturally organizes the decay rates of density fluctua- 
tions into short time scales at large wavevectors and long 
time scales at small wavevectors. As we comment fur- 
ther below however, there are technical obstacles to be 
overcome before this RG can be carried out. 

We take the disorder to be gaussian distributed, but 
not necessarily white noise, i.e. possibly smoothed. We 
then have 



««(q)» = o 



and 



ll l \ l 2tTQ! 2 V 2 _lfV( a 2 -11 r 

<(v(q)w(q')» = L 2 e 2?l < 



which in real space translates into 
«V(r)V(rO» = 



2 2 

a v 



P(a 2 



1 



■ exp 



q+q',o 5 



-'12 



2£ 2 (a 2 - 1) 



(23) 



(24) 



(25) 



Here v denotes the strength of the disorder potential and 
a is a dimensionless smoothness parameter. In the limit 
of a distribution which is extremely smooth (a — > oo), 
the one-particle electronic density of states approaches a 
gaussiano 



Pa=oo(E) = 



1 



(2tt) 3 / 2 £ 2 v 



exp 



2> 



(26) 



An integration over all energies E gives the number of 
states in the LLL divided by the sample area N/L 2 which 
is l/(27r^ 2 ). For a = 1 the disorder distribution goes over 
to the uncorrelated white noise distribution_£pr which 
Wegner has determined the density of states.c3 At E = 
E c it is given by 



p a =i{E c 



V2 
ir 2 £ 2 z 



(27) 



III. SELF-CONSISTENT BORN 
APPROXIMATION 

We next calculate Tl(q, w) Eq. ( p2| ) in the self-consistent 
Born approximation. We define the complex self-energy 
S(g, u>) = Tin(q,oj) + i'Ei(q,u>) for the propagator 



by setting 



tl(q,w) = ne 2 e^ 2q2 ii(q,Lu) 



U(q,u) 



uj + id — S(g, u>) 



(28) 



(29) 



Within the self-consistent Born approximation the self- 
energy is given by the expression 

s s ( g , w ) = ^<(/: q , q+p /: q+p , q ))n B (|q + p|, w ) . (30) 



In contrast to standard many-body perturbation theory 
the right hand side of this expression does not contain an 
energy sum. In this approximation all non-crossing dia- 
grams for the propagator U(q, lo) are summed, as shown 
111 Fig. § In this fi gure a thick solid line stands for 
tl B (q,w) and a thin solid line indicates the "noninteract- 
ing" propagator TL°(q,uj), which is given by Eq. (p9|) with 
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E(<7, lo) = 0. ft is independent of q. The consequences 
of this for a perturbative treatment will be discussed in 
the next section. In Fig. [I] the vertex with an incom- 
ing and an outgoing solid line and a dashed line stands 
for a matrix element £ qq ' of the Liouvillian. The disor- 
der average introduces "contractions", i.e. connections, 
between the dashed lines. In general the Hartree terms 
are included in the partial sum, but as indicated in Fig. 
[I], they vanish because of the q A p term in the matrix 
elements of the Liouvillian. 

Using the distribution introduced in the last section 
[see Eq. (p4|)] and the definition of the Liouvillian matrix 
Eq. (pit) we obtain the self-consistency equation 



T B (q,co) 



2ira 2 v 2 



Ii-Pl 2 4sm 2 (f qAp) 



p 



+ iS -T b ( p ,lo) 



(31) 



The strength of the disorder v can be scaled out of this 
equation by replacing S s — > KT B jv and lo — > hoj/v. 




U 2 Imn B (?,w) 



Sf(9,«) 



[lo - Eg (q,u) 



(32) 



Eq. (|3l]) can be solved numerically by iteration. Follow- 
ing Eq. (|l0|) the best way to extract the diffusive prop- 
erties is a "scaling" plot in which (for fixed v and a) 
— Tilo £ 2 Im H B (q, lo) is plotted as a function oiqL^ for dif- 
ferent small q and lo. Such an evaluation shows that on 
this level of approximation — TiLot 2 IvnXl B (q, lo) is a func- 
tion of qLu only and thus does not display a sign of 
the prefactor (q^o) x l v discussed in connection with Eq. 
(|l0|). Furthermore ImII B ((7,w) only shows normal diffu- 
sion with a diffusion constant Do which for q — > and 
lo — ► is independent of qL u . We thus conclude that 
(as expected) the occurrence of the critical exponents v 
and r\ is a higher order fluctuation effect. For lo — > 0, 



Y, B (q,Lo) goes to zero for all q. Thus D is given by 



Do 



lim lim T? (o, lo) I q 2 



(33) 



Because of the scaling property discussed following Eq. 
( pl| ) D is proportional to w. As shown in Fig. || D 
also depends on the smoothness a of the disorder. Be- 
tween a = 1 (white noise) and a = 2, Do changes 
by approximately ten percent. For a > 2 the a de- 
pendence is extremely weak and for a — > oo, Do sat- 
0.828 i^ 2 /fi. For a = 1 we find 

0.965 vi 2 /h. 

Using the Einstein relation for the conductivity and 
Eqs. (piq ) and (|27|) we obtain 



urates at Dq =cc 



+ 



+ 



— + 



FIG. 1. Partial sum of all non-crossing diagrams of the 
propagator Tl(q,Lo). For details see the text. 



0.330- 



(34) 



and 



1.00 



0.95 



a 0.90 



D 



0.85 



0.80 



FIG. 2. Diffusion constant Do as a function of the smooth- 
ness of the disorder a. 

As explained in the introduction the diffusive proper- 
ties can be read off from the small q and lo limit of the 
imaginary part of n. For q — > we have tl(q, u>) — > 
h~£ 2 H(q, co) and can thus write 



T^" 1 « 0.869 — 
xx h 



(35) 



If one is interested in the large a limit it might be tempt- 
ing to expand the sine in Eq. (^TJ) , as only small p con- 
tribute to the sum due to the exponential function. An- 
ticipating that for small q the self-energy is quadratic in 
q the ansatz T, B (q, lo) = —iq 2 D seems to be plausible. 
Then the self-consistency equation can be solved analyt- 
ically leading to D = (l/V2)v£ 2 /h w 0.707v£ 2 /h. A 
comparison with Dq = °° discussed above shows that this 
procedure does not give the correct large a value for Do. 
This is due to the fact that in the exact solution of Eq. 
( pl| ) the range of q values over which S s ((7, lo) can be ap- 
proximated by a purely quadratic function in q shrinks as 
Xjct. Thus in the limit a — > oo it would be necessary to 
include higher order terms in the expansion of TP (q, lo) 
in order to reproduce the numerical result in Eq. (|34). 
Note that a"^ 1 obtained above is independent of the 
correlation parameter a as it should in the limit a> 1. 
Since the exact conductivity is universal, the present re- 
sult is a considerable improvement over the traditional 
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self-consistent Born approximation result fpi which the 
conductivity vanishes like a -1 in this limit E3 

In a previous numerical studyliU it was found that 
(T xx = (0.54 ± 0.04)e 2 /ft. independent of the smoothness 
of the disorder. The results for cr xx obtained within our 
approach are of the same order of magnitude as the one 
calculated using purely numerical methodstll but in con- 
trast to this one our results depend on a. This is due 
to the fact that we have calculated Do within the self- 
consistent Born approximation but included in the Ein- 
stein relation the exact density of states at the critical 
energy. 

Using our approach of calculating the disorder aver- 
aged one-particle correlation function for the dynamical 
system described by the Liouvillian, we observe normal 
diffusion already at the level of the self-consistent Born 
approximation. In the usual fermionic picture of non- 
interacting electrons in the presence of disorder and a 
magnetic field, much more elaborate techniques, as e.g. 
Borel resummation, instanton methods, the replica trick, 
and the supersymmetry method, are used to obtain sim- 
ilar results. IH In particular, in the more traditional ap- 
proaches, diffusion is not obtained at the saddle point 
level and it is necessary to include gaussian fluctuations 
(i.e. sum ladder diagrams) to obtain diffusion. Because 
we deal directly with the density itself, we obtain diffu- 
sion even at the saddle point level. 



IV. FIELD THEORETICAL APPROACH 

To go beyond the self-consistent Born approximation 
it might prove advantageous to bring our approach into 
a field theoretical framework. This is what we will do in 
this section. In reformulating n(g, u>) using field theoret- 
ical methods we use the gaussian integral identity 



1 



to + iS — C 



where 



and 



Oq^q) 



Stf, = -i 



1 

Z 



^ [u + iS- £] k)k / -0k' 

k,k' 



(36) 



(37) 



(38) 



The V'q denote complex (bosonic) fields and Z is given 
by 



z 



VipVtp e" 6 "* 



(39) 



In order to ensemble average over the disorder we intro- 
duce additionaLGrassmann variables to represent 1/Z as 
a path integraltl 



where 



S{ = -i}; 6c [w + i6- £]k,k' £k 

k,k' 



(40) 



(41) 



One can then carry out the ensemble average over the 
gaussian distributed disorder and obtains the generalized 
functional 



(42) 



where 



F(uj) = ^ [i~ iuJ + <5k) V'kV'k + (-i<*> + <5k) £k£k 



+ X! «£kk'£pp')) [ipkfppfpp'^k' 

k.k' p,p' 

+2M>ut P t v > ■ (43) 

Here we have let 5 — > 8^ so that we can generate the 
correlation functions by 



(Cq£ q ) w = <0q^q) w = ~ 



dZ(oj) _ dZ(u) 



d5 



(44) 



Once the disorder averaging is done we finally obtain 

F(u>) =~i^2 [(^ + ^q)^qV>q + ( w + *'<5q)£q£q 

q 

+ J! /C 1 ' 2 ' 3 ' 4 ) [V'qiV'qaV'qsV^ 
q1.q2.q3.q4 

+20 qi ^q 4 Cq 2 Cq 3 +^^^3^4] > ( 45 ) 

with 



/(l, 2,3,4) = - 



TroV i^ a2|qi _ qJ 2 , . ft 



4 sin ( —qi A g 4 



x sin [ —q 2 A q 3 ) ^ qi+q2iq3+qi . (46) 

In contrast to standard many-body theory the action 
Eq. ( p5| ) does not contain a sum over the frequency, to 
only enters this equation as an external parameter. As 
already discussed in the last section the noninteracting 
propagator (v — 0) is given by (to + id)^ 1 and does not 
depend on q. Thus a perturbation theory or RG pro- 
cedure can only be set up after a q dependent propaga- 
tor has been generated by self-consistently summing up 
an entire class of diagrams, as e.g. the non-crossing dia- 
grams in Sec. III. Furthermore the interaction / in Eq. 
( ^6| ) has an unusual momentum dependence compared 
to standard standard </> 4 theory of critical phenomena: It 
vanishes if one of the goes to zero and is periodic in 
the momenta. 



G 



Using the field theoretical approach we can reproduce 
the approximation discussed in Sec. Ill, which is usually 
called self-consistent mean-field or saddle point approx- 
imation in the present context. In the absence of sym- 
metry breaking, the middle of the three quartic terms in 
the action cannot contribute to the saddle point solution 
since its coefficient vanishes for qi = q4 and q2 = q3. 
Hence we can deal separately with the bosonic and the 
fermionic variables when discussing the saddle point so- 
lution. By performing the usual pairing of the fields in 
the quartic interaction term at the mean-field level we 
have 

+in MF (g 2 ,c i ;)5 q2 , q4 ^ qi Vq 3 . (47) 

Thus we can write 

F MF H = ]T ^ [-ioj + S + z£ MF (g, W )] ^q , (48) 



and use this in calculating 

zn MF ( g ,c)EE<^X F = 



ifl B (q,tu), (49) 



LU + i5-Z MF (q,uj) 

which reproduces the self-consistency Eq. (|l]) for the 
self-energy. 

At present we do not know how to evaluate the corre- 
lation function beyond the self-consistent mean-field ap- 
proximation in a controlled way. However, we hope that 
in the future it will be possible to analytically extend our 
results. 



V. NUMERICAL RESULTS 

In this section we will numerically calculate ImII(g, to) 
by exact diagonalization and verify the scaling hypothesis 
stated in Eq. (10). We closely follow th&. procedure and 
notation used by Chalker and Daniell.EHl Motivated by 
Eq. (B) we define 



S(r,u) = -— Imn(r,w) 



= (CY S(hco + E i - ^)^(0)^(r)^(r)V*(0))) . (50) 



The single particle wave functions ipi ( r ) can be expanded 
in the basis of the elliptical theta functions <^> m (r) 



N 

t/>i(r) = ^2 ai{m)(j) m (r) 

rn — 1 



(51) 



^ X ^ = ~rfTTT^ exp(iX m ^ s y/£ 2 ) 

xexp[-(x-X rn , s ) 2 /(2£ 2 )] , (52) 



and 



27T 2 

X m ,s — m — * + S L . 
L 



(53) 



Then the Fourier transform of Eq. ( pOj ) can be written as 

1 1|2J 



S(q,u) 



2n£ 2 N 2 



(54) 



hi 



where 



Qa(k,l) = N 



N 



x exp 



2J a,i(m)a*([m - I}) 

m=l 

2 



(55) 



and q = (2n/L)(k, I) = y -psjjik, I), with k, I integer. In 

Eq. ( |55| ) [m+Z] is defined as being m+l for 1 < m+l < N 
and m + I ± N otherwise such that 1 < \m + I ± N\ < 
N. In the numerical calculation we replace the delta 
function in Eq. ( |54| ) by a sharply_Deaked gaussian 8 1 {x) cx 
exp[— x 2 /j 2 ] with a broadeningEj 7 = 0.64 u/iV which is 
of the order of the level spacing. We then have 



S{q,u) 



1 



2ir£ 2 N 2 



q K(q,u), 



(56) 



with 



K(q,u;) = ^-^ ^-U-. (57) 

This function is suitable for a numerical investiga- 
tion.!^ We restrict ourselves to a white noise disorder 
distribution (a = 1). We calculate K(q,io) for values of 
2 < k 2 + I 2 < 25 and Two = jn, with 3 < n < 23, where 
the limits have been chosen such that L^ 1 < q < I" 1 
and Tiuj -C v but Tluj greater than the level spacing of the 
finite size system. The system sizes range from N = 200 
to N — 2000, and the number of disorder realizations are 
500 or 100 depending on the system size. All values of 
K(q,uj) were determined to an accuracy better than 1% 
in the disorder averaging. 



where 



7 



1/3 

- * 1/4 



3 
cf 

£■ 10" 

£ 

s 

_l 

CT 
< 




10" 



10" 



n[1+(k 2 +l 2 ) 2 /n 2 ] 1/2 /N 



FIG. 3. %LoK{q,Lo) for different fixed (gl^) 2 oc (fc 2 + l 2 )/n 
and small g and a; as a function of n[l + (fc 2 + Z 2 ) 2 /n 2 ] (1/2) /AT 
on a log-log scale. For clear comparison each data set has 
been multiplied by a constant factor A(qL u} ) and is labeled 
by the ratio (fc 2 + l 2 )/n. 

For a fixed and small value of qL^ (so that we are 
in the range of normal diffusion) and q, to — > we expect 
from Eq. ( p^| ) that Tiu;K(q, u>) scales as (Hlu/vY^ oc (q£) ~ . 
The scaling hypothesis is illustrated in Fig. where we 
plot A(qL UJ )l%LoK{q 1 uo) for fixed ratios of (qL^) 2 oc (fc 2 + 
Z 2 )/n as a function of n^/l + (fc 2 + l 2 ) 2 /n 2 /N on a log- 
log scale. Here each curve is multiplied by a constant 
factor A(qL LU ) [different for each (fc 2 + l 2 )/n ratio] to 
make the comparison of the different lines easier. Also 
the factor y/T+ (fc 2 + I 2 ) 2 /n 2 multiplying n/N oc to is 
used such that the curves line up horizontally. The fact 
that data calculated for different system sizes fall onto 
the same curve [for a fixed ratio of (fc 2 +l 2 )/n)] indicates 
that the limits chosen above for fc, I, and n do avoid large 
finite size effects. On the log-log scale the different data 
sets fall onto straight lines and can be fitted by power- 
laws (solid lines in Fig. ||). 

The localization exponent v extracted from the slope of 



the lines in Fig. 
{qL^f in Fig. p 



I is shown as a function of (fc 2 + l 2 )/n oc 
Within our error bars and for the qL^ 



considered, v is a constant. Its value v — 2.33 ± 0.05 
is in excellent, agreement with previous finite-size scaling 
studiesElid'ErtHl and strongly supports the scaling hypoth- 
esis Eq. (|To|). The fact that the lowest (fc 2 + l 2 )/n points 
seem to be moving upwards in Fig. ^ is an indication that 
there are still some finite-size effects for the low values of 
(k 2 + l 2 ). In contrast to previous numerical studies we are 
able to obtain information about the critical exponent v 
from systems of finite size without doing finite-size scal- 
ing. 



3.0 



2.5 



> 2.0 



1.5 



1.0 



feH^ f I 1 I 



0.0 



0.1 



0.2 



0.3 



(k 2 +l 2 )/n 



FIG. 4. Localization exponent v obtained from Fig. 3. The 
solid line corresponds to the average v. 



VI. CONCLUSION 

We have presented a new analytical and numerical ap- 
proach to the localization-delocalization transition in the 
LLL of the IQHE. By using the closed Lie algebra of 
the density operators in the LLL we are able to write 
the equation of motion for the densities in a closed form 
which can be solved formally. Using the solution of the 
equation of motion for the projected densities we can 
express the integrated spectral function J dES(q,uj; E) 
as the disorder averaged density of states of a dynam- 
ical system with a novel action. We show analytically 
that the self-consistent mean-field approximation of the 
integrated spectral function yields normal diffusion but 
it misses the critical scaling. However, it is encouraging 
to note that even at this level of approximation the lon- 
gitudinal conductivity is inj-approximate agreement with 
previous numerical studies J13 Finally, using exact diago- 
nalization, we are able to extract the localization critical 
exponent v from the integrate d sp ectral function by us- 
ing the scaling hypothesis Eq. ( |10| ) , without having to do 
finite-size scaling. We obtain v = 2.33. ± 0.05 in excellent 
agreement with previous studies Jjlid We hope that in 
the future it will be possible to extend our approach be- 
yond the self-consistent mean-field level and analytically 
extract information about the critical exponent v. 
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